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Abstract 

Angular distribution of the 12 C( 7 Li, 6 He) 13 N reaction at E( 7 Li) = 44.0 MeV was 
measured at the HI-13 tandem accelerator of Beijing, China. Asymptotic normal- 
ization coefficient (ANC) of 13 N -> 12 C + p was derived to be 1.64 ± 0.11 fur 1 / 2 
through distorted wave Born approximation (DWBA) analysis. The ANC was then 
used to deduce the astrophysical S(E) factors and reaction rates for direct capture 
in 12 C(p, 7) 13 N at energies of astrophysical relevance. 
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measured cr(6), DWBA analysis, deduced asymptotic normalization coefficient, 
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1 Introduction 



In stellar evolution, the 12 C(p, 7) 13 N reaction plays a key role for the following 
reasons. Firstly, it is the first reaction in the Carbon-Nitrogen-Oxygen (CNO) 
cycle which dominates the energy production in stars with masses heavier than 
1.5 M Q [1,2]. Secondly, the 12 C(p, 7) 13 N(/3+) 13 C reactions can enhance the 13 C 
abundance [3], and thus influence the 12 C/ 13 C ratio which is thought to be an 
important measure for stellar evolution and nucleosynthesis [4]. Thirdly, the 
supply of 13 C by the 12 C(p, 7) 13 N(/? + ) 13 C reactions is also important for the 
13 C(a,n) 16 neutron source in the asymptotic giant branch (AGB) stars. A 
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recent calculation with parametric one-zone nucleosynthesis showed that the 
more 13 C supply is needed at the end of the CNO cycle in solar-metallicity 
stars [5]. In view of the above mentioned significance, it is highly desired to 
carefully investigate the 12 C(p, 7) 13 N cross section at energies below 1.0 MeV 
for the astrophysical interest. 

The 12 C(p, 7) 13 N reaction has been studied over a wide energy range down to 
about -E c .m. = 70 keV [6,7,8,9,10,11,12,13] since 1950. However, there exists 
obvious discrepancy in the experimental data at the lower energies. In the 
even lower energy range, no experimental data is available so far, and the 
astrophysical S(E) factors can only be derived by the extrapolation which 
results in uncertainty inevitably. In the energy range of E c . m . < 200 keV, the 
12 C(p, 7) 13 N reaction is dominated by the tail of the s-wave capture into the 
broad l/2 + resonance at E T = 421 keV. Although the contribution from direct 
capture is believed to be much smaller than that from the resonance tail, the 
interference between the two processes can lead to a considerable variation of 
S(E) factors, since both of them proceed via s-wave and then decay by El 
transitions. As a result, the S(E) factors either increase by the constructive 
interference or decrease by the destructive one. Thus the reliable experimental 
data on the direct capture are needed to derive the S(E) factors, particularly 
at energies of -E c . m . < 200 keV. A practicable scheme to deduce the S(E) 
factors for the direct capture in 12 C(p, 7) 13 N is combining the asymptotic 
normalization coefficient (ANC) of 13 N — > 12 C + p and R-matrix approach 
[14,15], the ANC can be deduced from the angular distribution of one proton 
transfer reactions. The ( 7 Li, 6 He) reactions are considered to be a valuable 
spectroscopic tool because the shapes of their angular distributions can be 
well reproduced by the distorted wave Born approximation (DWBA) [16]. 
Thus, the 12 C( 7 Li, 6 He) 13 N reaction is used to extract the nuclear ANC of 13 N 
— > 12 C + p in the present work. 

The 12 C( 7 Li, 6 He) 13 N angular distribution was measured at E( 7 Li) = 44.0 
MeV. The spectroscopic factor and ANC were derived based on DWBA anal- 
ysis, and then used to calculate the astrophysical S(E) factors and rates of 
12 C(p, 7) 13 N direct capture reaction at energies of astrophysical interest with 
the R-matrix approach. We have also computed the contribution from the 
resonant capture and the interference effect between resonant and direct cap- 
tures. 



2 Measurement of the 12 C( 7 Li, 6 He) 13 N angular distribution 

The experiment was performed at the HI-13 tandem accelerator of Beijing, 
China. A carbon target in thickness of 39.0 /xg/cm 2 was bombarded with 
the 44.0 MeV 7 Li beam in intensity of about 100 pnA. The 6 He ions from the 
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Fig. 1. Angular distribution of the 12 C( 7 Li, 6 He) 13 N reaction at E('Li) = 44.0 MeV. 

12 C( 7 Li, 6 He) 13 N reaction were analyzed with the Q3D magnetic spectrograph. 
In order to attain a good angular resolution, the solid angle for the reaction 
products was set to be 0.23 msr. A two dimensional position sensitive silicon 
detector (PSSD) was placed at the focal plane of the spectrograph, which 
assured the full detection of 6 He ions emitted within the solid angle. During the 
measurement, a Faraday cup placed at 0° in the reaction chamber was utilized 
to record the beam current which served as the normalization standard in 
determining the absolute cross section, while two independent silicon detectors 
(SSDs) at ± 30° were used for both monitoring the beam balance and relative 
normalizing the measured cross sections. The normalization was checked for 
several angles from 23° to 40° with the elastic scattering data of 36.0 MeV 7 Li 



The 12 C( 7 Li, 6 He) 13 N differential cross sections were measured in the angular 
range of 7° < 6^ < 23°, corresponding to 11° < 9 c ni < 37°. The measured 
angular distribution is shown in Fig. 1. The experimental errors are from the 
uncertainties of statistics, target thickness (5%) and solid angle (6%). 



3 Determination of the 13 N Nuclear ANC 

The spins and parities of 12 C and 13 N ground states are + and l/2~, re- 
spectively. The 12 C( 7 Li, 6 He) 13 N cross section is dominated by the (0 + ,0) — > 
(l/2~,l/2) transition, thus only lpi/ 2 orbit in 13 N can be populated. If the 
reaction is peripheral, the differential cross section can be expressed as 



on 



C [17]. 




(1) 
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Table 1 

Optical potential parameters used in the DWBA calculations, where U, W are in 
MeV, r and a in fm. 



Channel 


U v 




aR 


W v 


ri 


ai 


r c 


7 Li + 12 C 


194.55 


0.50 


0.82 


7.66 


1.31 


0.75 


1.30 


6 He + 13 N 


170.50 


0.79 


0.67 


10.70 


1.31 


0.75 


1.30 



where (^) eX p and (^)dwba denote the measured and calculated differential 
cross sections respectively. C™y 2 and Ci Li are the nuclear ANCs of 13 N — > 12 C 
+ p and 7 Li — > 6 He + p, b^y 2 and 6? Li being the single particle ANCs of the 
bound state protons in 13 N and 7 Li, which can be calculated with the bound 
state single particle wave function and Whittaker function at larger radius. 
The ratio of C/b is so called spectroscopic factor. The proton spectroscopic 
factor of 7 Li was derived to be 0.41 ± 0.05 in our previous work [18]. The 
nuclear ANC of 13 N ground state can then be extracted by normalizing the 
theoretical differential cross sections to the experimental data via Eq. (1). 

The angular distribution of the 12 C( 7 Li, 6 He) 13 N reaction was calculated with 
the DWBA code FRESCO [19]. The bound state wave function was obtained 
by solving the Schrodinger equation using a Woods-Saxon potential with stan- 
dard geometrical parameters (ro = 1.25 fm and a = 0.65 fm), the potential 
depth was adjusted so as to reproduce the observed binding energy of the va- 
lence proton. The optical potential parameters for both the entrance and exit 
channels were extracted by fitting the angular distributions of elastic scat- 
tering with the code SFRESCO [19]. The 7 Li + 12 C elastic scattering was 
measured in the present work. Since there is no experimental data for 6 He 
elastic scattering on 13 N , the angular distribution of 6 Li + 13 C elastic scatter- 
ing at E( 6 Li) = 28.0 MeV [20] was used to extract the potential parameters 
of the exit channel. The extracted potentials are listed in Table I. 

Generally, the spectroscopic factor or nuclear ANC is determined by fitting 
the theoretical calculations to the experimental data at the first peak in the 
angular distribution for the forward angles [21], since the experimental differ- 
ential cross sections for the backward angles are very sensitive to the inelastic 
coupling effects and other high-order ones, which can not be well described 
theoretically. In the DWBA calculation, the differential cross sections at three 
forward angles were used to extract the ANC. The normalized angular dis- 
tribution is also presented in Fig. 1. The nuclear ANC and the spectroscopic 
factor for 13 N -> 12 C + p are deduced to be 1.64 ± 0.11 fm^ 1 / 2 and 0.64 ± 0.09, 
respectively. The spectroscopic factors extracted by different experiments and 
the theoretical values are listed in Tab. 2. The spectroscopic factor obtained 
in our work agrees with the theoretical ones reported in Refs. [22,23] and the 
experimental results given in Refs. [26,28,29,37]. The Nuclear ANC of 13 N — > 
12 C + p from this work is in good agreement with the value of 1.65 ± 0.20 
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Table 2 

The theoretical and experimental proton spectroscopic factors in the 13 N ground 
state. 





Experiments or theory 


Year 


Reference 


0.61 


theory 


1967 


[22] 


0.56 


theory 


1969 


[23] 


0.78 - 1.35 


l2 C(d,n) 


1970 


[24] 


0.74 


12 C(d,n) 


1971 


[25] 


0.53 ± 0.12 


12 C(d,n) 


1972 


[26] 


0.70 - 1.48 


12 C( 3 He, d) 


1969 


[27] 


0.56 - 0.78 


12 C( 3 He, d) 


1976 


[28] 


0.68 ± 0.12 


12 C( 3 He, d) 


1976 


[29] 


0.81 ± 0.12 


12 C( 3 He, d) 


1979 


[30] 


0.48 ± 0.12 


12 C( 3 He, d) 


1980 


[31] 


1.34 


12 C{a,t) 


1969 


[32] 


0.91 


12 C(a,t) 


1972 


[33] 


0.72 


12 C( 7 Li, 6 He) 


1979 


[34] 


0.38 ± 0.05 


12 C( 7 Li, 6 He) 


1986 


[35] 


0.64 ± 0.09 


12 C( 7 Li, 6 He) 


2008 


present work 


0.25, 0.40 


12 C( 10 B, 9 Be) 


1974 


[36] 


0.62 


12 C( 14 N, 13 C) 


1975 


[37] 


0.29, 0.40 


12 C (16 0) 15 N ) 


1979 


[38] 



fm- 1 / 2 extracted from the 12 C( 10 B, 9 Be) 13 N reaction [39]. 



4 Astrophysical S(E) factors of the 12 C(p, 7) 13 N reaction 

Following the approach used in our previous work [15], the astrophysical S(E) 
factors of the 12 C(p, 7) 13 N reaction were calculated by the R-matrix method. 
For the radiative capture reaction B + b — > A + 7, the cross section to the 
state of nucleus A with the spin Jf can be written as [14,15] 

a Jf = J2 a JiJ f i ( 2 ) 
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where J, denotes the total angular momentum of the colliding nuclei B and b 
in the initial state, Jb and Jb are the spins of nuclei b and B, and I, k and 
Zj are their channel spin, wave number and orbital angular momentum in the 
initial state, respectively. Un i j s j l is the transition amplitude from the initial 
continuum state (Jj,/,/j) to the final bound state (J/,/). In the single-level, 
single-channel approximation, the resonant amplitude for the capture into the 
resonance with energy E Rn and spin Jj, and subsequent decay into the bound 
state with the spin J/ can be expressed as 

tjR - ? - r ^,-^,) fe(^ r fe( £; )] 1/2 m 
""* J <~ E-E Rn+ iTjj2 ■ (4) 



Here it is assumed that the boundary parameter is equal to the shift function 
at resonance energy, and 4>u is the hard-sphere phase shift in the /jth partial 
wave, 

rFj.(fc,r c )i /r \ 



where F 2 and Gf. are the regular and irregular Coulomb functions, r c is the 
channel radius. The Coulomb phase factor uj Xi is given by 

k 

u k = arctan(— ), (6) 

n=l U 



where rji is the Sommerfeld parameter. r b } z . (E) is the observable partial width 
of the resonance in the channel B + b, F^j (E) is the observable radiative 
width for the decay of the given resonance into the bound state with the spin 
Jf, and Fj i pa J2 ^lii ls the observable total width of the resonance level. The 

energy dependence of the partial widths is determined by 

r &(^) = ^Pr r &(^J (7) 



and 



rf Jf (E) = (^^) 2L+1 ^ Jf (E Rn ), (8) 
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where T^j^^E^) and T^ i j f (Eji n ) are the experimental partial and radiative 
widths, Ef is the proton binding energy of the bound state in nucleus A, and 
L is the multipolarity of the gamma transition. The penetrability Pi { (E) is 
expressed as 

Pk{E) = m,r e ) + GUk,r c y (9) 



The nonresonant amplitude can be calculated by 



U? h % Jt = -(2) 3 /V^+VK-^^ + (-1)^1 ( Vc) L+1/2 



x - 



{L + 1] l L + ^ {2L I ^ C JfIlf F h (k,r c )G h (k,r c 



xW^Kr^^ikOLOllfOpiLlfJJ^^J'^klf), (10) 

where fc 7 = (E + £f)/hc is the wave number of the emitted photon. Cj f n f is 
the nuclear ANC of 13 N — > 12 C + p, and // are the wave number and rela- 
tive orbital angular momentum of the bound state. Wi(2nr) is the Whittaker 
hypergeometric function with k = y/'2fMBb£f- (h0L0\lf0) and U(LlfJiI;liJf) 
are the Clebsch-Gordan and Racha coefficients, respectively. Jiiklf) is the 
integral expression defined as 



oo 



j' (Hf) - _^ / dr ^ W^ r F^r) _ G h {k,r) ] 



The non-resonant amplitude contains the radial integral ranging only from 
the channel radius r c to infinity since the internal contribution is contained 
within the resonant part. Furthermore, the R-matrix boundary condition at 
the channel radius r c implies that the scattering of particles in the initial state 
is given by the hard sphere phase. Hence, the problems related to the interior 
contribution and the choice of incident channel optical parameters do not 
occur. Therefore, the direct capture cross section only depends on the ANC 
and the channel radius r c . Using the experimental ANC (1.64 ± 0.11 fm~ 1//2 ) 
from the present work, the non-resonant 12 C(p, 7) 13 N cross sections vs. E c . m . 
were calculated, as shown in Fig 2. In the calculation, r c was taken to be 5.0 
fm following the previous works in Ref. [15] and [40]. 

The astrophysical S-factor is related to the cross section by 

S(E) = Ea(E) exp{E G /Ey/ 2 , (12) 
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Fig. 2. The non-resonant cross sections of the 12 C(p, 7) 13 N reaction computed with 
the ANC derived from the present experiment. 
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Fig. 3. Astrophysical S(E) factors as a function of -E c . m . for the 12 C(p, 7) 13 N reac- 
tion. The dashed line and the solid line are the contributions from the direct proton 
capture and the total S(E) factors, respectively. The experimental data are taken 
from Refs. [6,7,8,9,12] 

where the Gamow energy Eq = 0.978ZfZ%fi MeV, /i is the reduced mass of the 
system. Using the resonance parameters (E R =421 keV, T tot (E R ) = 36.5 keV, 
and T^(E R ) = 0.67 eV) from Ref. [41], the S-factors for direct and resonant 
captures can be then derived. 

Since the incoming angular momentum (s-wave) and the multipolarity {El) 
for the direct and resonant capture 7-radiation are identical, there is an in- 
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terference between the two capture processes. Following the Ref. [9], the total 
S-factor is calculated by 

Stot(E) = S dc (E) + S res (E) ± 2[S dc (E)S res (E)] 1/2 cos(5 r ), (13) 



where 5 r is the resonance phase shift, given by 

r P (£) 



arctan 



2{E - E r ) 



(14) 



The experimental results from the direct measurement of 12 C(p, 7) 13 N reaction 
[9] show that the interference between the resonant and direct captures is 
constructive below the resonance energy, and destructive above it. Based on 
this interference pattern, the present total S-factors are then obtained, as 
shown in Fig. 3. One can see that the present S-factors of 12 C(p, 7) 13 N are in 
good agreement with the experimental data from Ref. [42] and the references 
therein. Very recently, N. Burtebaev et al. [12] measured the cross sections of 
the 12 C(p,7) 13 N reaction at beam energies E p = 354, 390, 460, 463, 565, 750, 
and 1061 keV. They obtained ANC( 13 N) = 1.72 far 1 / 2 , Tj = 0.65 ± 0.07 eV 
and = 35.0 ± 1.0 keV by fitting their experimental data with R-matrix 
approach. Their results are in good agreement with ours. 



5 Astrophysical 12 C(p, 7) 13 N reaction rates 



The astrophysical 12 C(p, 7) 13 N reaction rate is calculated with 



N A{av) = iV,(^_i_ j S{E)eM - - JL] rfB , (15, 



where v = y2E//j,. Na and ks are Avogadro and Boltzmann constants re- 
spectively. The updated reaction rates are shown in Fig. 4, together with the 
previous ones from NACRE and CF88 compilations. The present reaction 
rates are larger than that of NACRE and CF88 by about 20% in the low 
temperature range of T 9 < 0.5, and very close to their average value in the 
temperature range of 0.5 < T 9 < 10.0. 

The total reaction rates as a function of temperature obtained in our work 
are parameterized with an expression used in the astrophysical reaction rate 
library REACLIB [43], 
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Fig. 4. (Color online) The 12 C(p, 7) 13 N reaction rates from present work as well as 
those from NACRE and CF88 compilations. 



Table 3 

The fitting parameters of the 12 C(p, 7) 13 N reaction rates. 



parameters 



value 



parameters 



value 



a 

do 
0,4 

a 5 

«6 



1.1397E+01 
-2.9913E-02 
-1.1507E+01 
7.6988E+00 
-3.4841E+00 
3.8122E-01 
4.5641E-02 



bo 
bi 
b 2 
b 3 
bi 
h 

ba 



4.0200E+01 
2.1886E+00 
-1.5188E+02 
1.1973E+02 
-3.7889E+00 
1.2086E-01 
-8.1517E+01 



N A (av) = exp(a + a{T 9 1 + a 2 T 9 1/3 + a 3 T 9 1/3 + a 4 T 9 

+ a 5 Tg /3 + a 6 InTg) + exp(6 + hTg 1 + b 2 T 9 1/3 
+ hT 9 1/3 + b 4 T 9 + 6 5 T 9 5/3 + b 6 lnT 9 ). 



(16) 



The value of fit parameters a _ 6 and &o-6 are listed in Tab. 3, and the fitting 
errors are less than 1% in the temperature range of 0.01 < T 9 < 10. 
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6 Conclusion and discussion 

The 12 C(p, 7) 13 N reaction is of considerably astrophysical interest. The cross 
sections in the energy range of astrophysical relevance is very small and dif- 
ficult to be measured directly. In the energy range of E cul _ < 200 keV, the 
12 C(p, 7) 13 N cross section depends on the tail of the s-wave capture into the 
broad l/2 + resonance at E T = 421 keV, the direct capture to the ground state 
and their interference. The determination of the resonant parameters for the 
l/2 + state and the ANC of 13 N ground state is helpful in extrapolating the ex- 
perimental data at high energies down to the energies of astrophysical interest 
(around 25 keV). 

We measured the angular distribution of the 12 C( 7 Li, 6 He) 13 N reaction at 
= 44.0 MeV, and deduced the nuclear ANC and spectroscopic factor 
for the 13 N ground state. The astrophysical S-factors and reaction rates of 
12 C(p, 7) 13 N are then extracted with the R-matrix approach. The S-factor at 
25.0 keV is found to be 1.87 ± 0.13 keV-b. The result is consistent with that 
of 1.75 ± 0.22 keV-b obtained by Burteaev et al. [12], and slightly higher than 
the earlier values of 1.45 ± 0.20 keV-b, 1.54 ± 0.08 keV-b and 1.33 ± 0.15 
keV-b reported in Refs. [9,10,11]. In order to clarify the deviation, the further 
study of the resonant parameters of the l/2 + state is needed. 
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